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Comparison of tho Nimbus-4 BUV Ozone Data 
Vith tha Ames Two-Dimensional Modal 

By U. J. Boruckl 1 ) and I. J. Ibarataln 2 ) 

Abstract — A comparlaon la made of the flrat two yaara of Nimbus-4 
backacatterad ultraviolet (BUV) otonc aeaaureaenta vith tho pradictiona of 
the Ames two-dimensional model. The ozone obaervatione used in thia atudy 
consist of the mixing ratio on the 1-, 2-, 5-, and 10-nb prasaure aurfaces. 
These data are zone and time averaged to obtain seasonal moans for 1970 and 
1971 and are found to show strong and repeatable meridional and aaaaonal 
dependencies. The model used for comparison with the observations extends 
from 80* N to 80* S latitude and from altitudes of 0 to 60 km with 5* hori- 
zontal grid spacing and 2.5-km vertical grid spacing. The chemical reaction 
and photolysis rate constants used In the model are those recommended In the 
report of the NASA Panel for Data Evaluation (1979). Chemical reaction 
and photolysis rates are dlurnally averaged, and the photodlsaociatlon rates 
are corrected for the effects of scattering. 

It is found that the large altitude, latitude, and aaaaonal changes In 
the ozone data agree well with the model predictions. Also shown are model 
predictions of the sensitivity of tha comparisons to changes In the assumed 
mixing ratios of water vapor, odd nitrogen, and odd chlorine, as well as 
to changes in the ambient temperature and transport parameters. 
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1, Introduction 

Because oione effects the aolar UV radiation that reaches the earth's 
surface and because osone influences stratospheric circulation patterns, it 
is important to determine the global distribution of osone and to understand 
the mechanisms that control the osone distribution. 

Model predictions of the variation of the osone distribution with 
altitude, latitude, and season allow us to comprehend the mechanisms that 
control the production, loss, and movement of osone. Hence, it is important 
to continuously compare model predictions with the measured distributions 
over as vide a range of conditions as possible to assess the depth of our 
understanding of these mechanisms and the validity of our predictions of the 
effects of anthropomorphic disturbances. 

As so often happens, the easiest measurements to make have been precisely 
the one that are the most difficult to SK>del, l.e., measurements at one point 
in space and time. However, with the advent of satellite measurements of 
osone by HEATH et al. (1973) and PRABHAKARA et si. (1973) , measurements of 
ozone distributions became available that are much more representative of 
seasonal and meridional means and are thus more amenable to comparisons with 
currently available models. These measurements are not yet completely satis- 
factory in that the BUV (backscattered ultraviolet) measurements do not cover 
the winter pole and the IRIS measurements do not give the vertical distri- 
butions of ozone. Further, these measurements are not entirely consistent 
they suffer from problems associated with mathematical difficulties in 
deconvolving the measurements to obtain the vertical dlatrlbutlona. Finally, 
the satellite data are available only for fairly short periods of time and 
hence might not be representative of long-term means. Nevertheless, the 
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satellite data offer a significant improvement over previous data in that 
their Internal consistency and large geographical coverage allow a more real* 
is tic determination of the vertical* meridional, and aaasonal distributions 
of osone for comparison with model predictions. 

It is the purpose of this paper to compare the Ames two-dimensional 
model predictions of osone concentrations with the BUV observations from the 
Nimbus-4 satellite and to Interpret the comparisons. Such a comparison 
should be viewed as an early attempt at model validation. As satellite data 
on other minor constituents (such as oxides of nitrogen, hydroxyl, and 
halogens) become available, it will become possible to make more complete, 
and thus more sophisticated, comparisons. 

2. Nimbus IV BUV Experiment 

The BUV experiment is on the Nimbus-4 satellite, which was launched into 
a circular polar orbit at an altitude of 1100 km. The 10* retrograde orbit 
ia sun synchronous, and every 107 minutes the three-axls-stablllsed satellite 
passes through the ascending node near local noon. The subsatellite point 
crosses the equator in Increments of 27* in longitude between successive 
passes and reaches a maximum latitude of 80*. The apatlal resolution of the 
Instrument is 230 km, or approximately 2* latitude (HEATH et al. , 1973). 

The Instrument la a double (tandem) Ibert-Fastle spectrophotometer in 
conjunction with a narrow band Interference filter photometer (HEATH et al.. 
1973). The spectrometer is designed to measure spectral intensities at 
twelve wavelengths between 2S5.5 nm and 339.8 nm with a 1-nm bandpass. The 
resulting light Intensity is measured with a hlgh-quantum-ef f lclency photo- 
multiplier tube. The absolute error of the intensity measurement is estimated 
at ilOX while the relative error is only 2X to 5X (FLEIG, 1980; HEATH et al. . 
1979). 
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The principle of the BUV experiment roots on the fact that ataoapharic 
radiance at a particular wavelength orlginatea in a moderately well-defined 
effective scattering layer. An outline of the principle of the experiment can 
be deduced froa Fig. 1. The Intensity of light reaching the observing instru- 
ment is proportional to the aolecular density In the scattering layer. By 
using the hydrostatic aquation and the Ideal gas law It is readily shown 
that absorption by oxone along the path shown in Fig. 1 is proportional to 
the pressure at the scattering level. 

Thus, we have the light reaching the receiver given by the following 
relation of proportionality 


I - .(l - T 0 ! ° x p)C8 x P> (1) 

where - ozone absorption coefficient 

Y^ • eeau ozone nixing ratio along the optical path 
0^ - scattering cross section 

P - pressure 

The paraaeter B is a constant of proportionality. It nay seen froa 
equation (1) that I goes to zero at the top of the atmosphere, where P - 0, 
and at the "bottan" defined by asxlnua penetration depth, wher 1 - Y Q ^oP - 0. 
Between those extreaes, I Bust go through a aaxlaal value, whose position nay 
be found by differentiating equation (1) with raspact to P and setting the 
derivative equal to zero. The pressura corresponding to naxlaun intensity 
is thus given by 


box 


*Vx 


( 2 ) 


It may be seen froa equation (2) that tha height of the maximum of the con- 
tribution function depends on the ozone absorption coefficient at the 
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wavelength X end on the Mount of otone preaent. A eet of repreeentetlve 
weighting function* la ahawn in Fig. 2. 

2. Bl fV Data Raduotim 


Following TWOMEY (1961) and MATEER (1972) we nap write the expreaelon 
for the backacatter radiance aa followa: 


X(X,e) 

s<^»Xp) 



♦ co* 2 e)s(v 4 
aec 8) (o^Xp 


V 


♦ B x P)|dP 


(3a) 

(3b) 


where I ■ backacattered radiance in the eatelllte nadir direction 
6 ■ eolar senlth angle 
X - wavelength of radiation 
v * the wave nunber, 1A 
Fq(X) • extraterrestrial eolar irradlance 
0 X ■ atnoepherlc acattering coefficient 
* oaone abeorptlon coefficient 
X ■ ozone eolunn denalty above praaeure P 


In equation (3) a plane parallel atnoephsre haa been aeaunad. Thua, addi- 
tional correction tome are needed for high aolar zenith angle*. 

In the upper stratosphere, abeorptlon by ozone la ouch atronger than 
Rayleigh acattering, and thua the acattering tern g^P nay be neglected 

conpared with the abaorptlon tern a.X . 

* P 
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equation (3b) 


uy thus bo rewritten no followot 


8(0.X p ) 



I«*p(-kX p )]dP 


(3c) 


vhoro tho prosauro Independent tana have bean gathered In the quantity 

k S (1 ♦ aec 6) 


Following TWOMEY (1961), the osona colunn density, X p , can be aade Into the 
variable of Integration. Thus the right aide of equation (3c) the Laplace 
transform of the derivative of P with respect to X p . Accusing the P(X p ) 
le bounded by an exponential less than exp(kX p ), equation (3c) now becomes 

S(v,X ) i Q(k) - kL[P(X )] (4) 

P P 

Since S(v,X ) Is a measured quantity, equation (4a) can be solved provided 
P 

that one can find a relationship between p end X . GREEN (1964) suggested 

P 

the following: 

X - CP l/j (5a) 

P 


P - C'°X 0 (5b) 

P 

where C and o are empirical parameters to be determined from the dsta. 

It is readily seen from equation (S) that a Is the ratio of ozone scale 
height to atmospheric scale height. This ratio la now asstned to be constant 
in the upper stratosphere. Substituting equation (3b) into aquation (4) and 
taking the Inverse transform (CHURCHILL, 1958, p. 324) yields the following 
expression: 

Q(k) - (Ck)“°r(o ♦ 1) (6) 
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The four shortest wavelengths of the experiment are now used to draw a 
straight line with elope -o on a log-lot plot, (lvlng osona overburden aa 
a function of ataospherlc praaaura. Dlf f erect iat log the raault glvaa osone 
nixing ratio. The above lnverelon nay be uaed for atnoapharle preaeurea leaa 
than 2 to 3 mb, l.e., altltudaa above 45 kn. For larger praaaura lave la, or 
lover altitudea, the aaeuaptlon that the loae of light lntanalty due to osona 
abaorptlon la much larger than that due to ecattering can no longer be nada, 
i.e. , the B^P tarn can no longar be naglactad conparod with Thus, 

it becoaes necesaary to solve aquation (3b) in its coop late form, which can- 
not be done by the slnple Laplace transform technique dlacusaad above. 

Mow define 


! (1 + aac 8}« x 

► 

6. 5 (1 ♦ aac 6)6. 

J A 

Then equation (3b) nay be written in the fora 

/ I -a X -8.P 
e 3 p e 3 dP 


(7a) 


(7b) 



The term SA\> At X ) has been written in aubacripted fora since it la not a 
J J P 

continuous function of frequency but la rather defined for a finite maber of 
frequencies. Equation (7b) la a fora of Fredhola's Integral aquation of 
the first kind (WHITTAKER and WATSON, 1958, p. 219), and it nay be written 

as 


♦(X) •» X 



b 

K(X,Od(OdC 


(7c) 
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FREDHOLM (1900) showed chat tha above Integral equation la the Halted fora 
(when 6 * 0) of the aqua t loot 


♦ (X) - Id r K(X.X)#(X > (8) 

P ^ P 9 9 


The above ay a ten of equations for 9(X ) (p"l» 2... 3) has a solution 

P 

If the detaralnant formed by the coefficients of #(X ) vanishes. The 

P 

expression for this determinant la given below froa WHITTAKER and WATSON (19S8). 


D(X) - 


1 - XdK(X l ,X l ) 
- X6K(X 2 ,X 1 ) 


- XdK(X x ,X 2 ) 
1 - X6K(X 2 ,X 2 ) 


- X6K(X l ,X n ) 

- XdK(X 2 ,Xn) 


( 9 ) 


X6K(X,X ) 
n i 


XdK (x n ,x 2 ) 


1 - XdK(X n ,X B ) 


Thus, the idealised problem reduces to finding eigen values of D(X). How- 
ever, the real situation Introduces special difficulties. Because of both 
noise in the experimental measurement and numerical Inaccuracies, equation (9) 
cannot be solved exactly. 

The Implicit Inversion problem now becomes one of minimising the quanti- 
ties: ’ 



( 10 ) 


where Is the measured radiance at the jth wave number and 1^ Is the 

radiance obtained by assuming a specific osone distribution. The quantity 

"j" varies over the number of wavelengths available for the experiment. In 

practice, 1 seldom exceeds 6 and Is sometimes less; J limits the 
' J sum J max 
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information that can be recovered in the oeone profile to that represented by 

a polynomial of order J . A further difficulty la cauaed by numerical 

Instabilities associated with the solution of equation (9). PHILLIPS (1962) 

has treated the case with numerical noise and developed a technique which is 

successful, but it imposes a smoothing constraint on the variation of X with 

P 

P in equation (9). This approach was further refined by TWOMEY (1963). 

In sunsary, the lnverslou procedure la as follovst 

1. An initial guess of X p (P) is constructed from a spline fit of the 
solution of the Laplace solution for altitudes above 2 mb, and a mean atmos- 
pheric profile is obtained from rocket data for altitudes below 2 mb. 

2. The backscatter equation (7b) is written in matrix form similar to 
equation (9) but incorporating the Phllllps-Twomey smoothing parameter, y. 

3. The numerical profile of X as a function of P is put into matrix 

P 

form, and the system is solved for the best profile consistent with experi- 
mental measurements. 

The detailed numerical procedure will be available In report form from 
FLEIC, et al. (1980). 

Analysis of the experiment indicates that the measurement error is 
between 22 and 52 (HEATH, et si., 1979). The inversion error ranges from 22 
at a 50- km altitude to 52 at a 30-km altitude (FLEIC, personal communication, 
1980). In addition, there is a sampling error of 1.52 to 2.52, depending on 
the amount of data that was accepted for the calculation of the seasonal 
means. Combining these errors gives a total experimental error of 32 to 7.52. 

A more complete analysis of the experimental errors can be obtained from the 
Goddard User's Manual for BUV Data by HUDSON et al. (1980). 
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4. Model dseeripticn 

The nodal la baaad on tha governing aquae Iona for tract conatltuanta 


3. ♦ 7 . ( V > ♦ ! • !, • s, 




( 11 ) 


( 12 ) 


where * concentration of the ith constituent 
• flue of the 1th conatltuent 

■ respective chealcal productlon^utd loaa rate tera 
H - ataosphcrlc acale height f* 

T ■ ataoapherlc acale teaperature 
e - vertical unit vector 

Z 

y - bulk velocity of the ataoephere 
K # • eddy dlffualvlty tenaor 

The chealcal rate equatlone are aolved by ualng an lapllclt technique. 
Oalttlng the transport ten froa equation (11) yields the finite difference 
fon: 


«s 


J+l i 

°1 " "l . s )♦» 

At 5 i 


( 13 ) 


which can be linearized by taking the first ten In a Taylor series expansion 
of about : 


n J+l _ J 

Oj * n 


At 


‘-.‘•Y&V' -<) 


( 14 ) 
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The aeabers of the set of aass conserve: Ion aqua t Iona ara coupled and raqulra 
aolutlon by aatrlx aathoda. Vlth tha lapllclc aethod, larga ciaa steps can 
ba takan wan chough tha aquation* ara stiff, 

f 

The tranaport coaputatlons ara tlae-spllt Into thoaa for vertical and 
horizontal advactlon and dlffualon. Solutlona for tha tranaport procaaaca 
arc obtained frots a nase-conservlng, forvard-tias, space-centered, finite 
difference fornulatlon of tha aquations. An implicit foniulation la used 

for tha horizontal diffusion and an explicit foraulation is used for all 

jb 

other tranaport processea. ' 

Spherical geoaetry is enployed in aquations (11) and (12) vlth the 

horizontal coordinate extending along a nerldian from 80* S to 80* N in 

5* intervals and the vertical coordinate extending frea the ground to 60 ka 

% 

in 2 . 5-lan Intervals. End boundaries ara takan at 80* S and 80* N because 
aeridlonal fluxes ara expected to ba snail at these latitudes. Tha and 
boundary conditions arc taken to be zero flux of all constituents across 
the vertical boundaries. The upper boundary conditions ate given by setting 
the vertical flux equal to sero for all species. Tha lover boundary condition 
for all apeclea except HN0 3 , H 2 0, M0 2 , Oj, MCI, H 2 C0, and H 2 0 2 is chaaical 
equilibrium Because HN0 3 , H 2 C0, MCI, and H 2 0 2 are water soluble, fixed nixing 
ratios are lapofife at the lover boundary that are consistent with the vash- 
out of these aolecules in the troposphere. The nun be r densities of N0 2 and 
N 2 0 are fixed at 1 « 10* cm” 3 and 7.5 * 10* 12 ca" 3 , respectively, at the 
lover boundary vhlle the nuaber density of 0 3 is fixed at 6 * 10 1 1 ca“ 3 . At 
each tine step, the nodel calculates the concentrations of 0 3 , 0( 3 P), 0( 1 D), 
N0 2 , NO, N0 3 , N 2 0, N 2 0 s, HN0 3 , HN0 2 , H 2 0 2i K0 2 , N, OH, H, Cl, CIO, MCI, C1N0 3 . 
CH 3 , HCO, H 2 CO. CH 3 0 2 . CH 3 0, CH„0 2 , do, CH 3 C1, CTC1 3 , and CF 2 C1 2 . The 
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profile* of N 2 . 0 2 , H 2 0, end H 2 are held fixed (at the experimentally 
determined valuta) during the calculation*. 

With » !n exceptiona t the value* of the photochemlcal-react ion-rate 
conatanta that ar* uacd are thoa* recomended in the report of the NASA Panel 
for Data Evaluation (1979) and are Hated in the Appendix. 

The chemical reaction and photolyaie rates were dlurnally averaged by 
using the techniques of TURCO and WHITTEN (1978) and COGLEY and BORUCKI (1976). 
The photodlsaociation rates have been corrected for the effect* of scattering 
by a technique based on the work of LUTHER and GELINAS (1976). 

The mean meridional circulation is obtained by the kinematic method from 
the steady-state equation of mass continuity in spherical coordinates: 

F^Te M < 55 CM 9 > * h <55) • 0 <15) 

where the overbar denotes a zonal and seasonal average 
p ■ density 

9 ■ meridional velocity component 
w • vertical velocity component 
R ■ Earth's radius 
6 » latitude 

z • altitude above mean sea level 

Equation (15) Implies the existence of a "stream function" for the total 
flux such that 


2*Rpv cos 8 • - 

(16a) 

2*Rpw cos 8 - ^ f 

(16b) 
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If the distributions of p and 9 are known, Y can be obtained from 
equation (16a) by vertical integration, and 0 can be obtained from equa- 
tion (16b). The 9 components were baaed on the obaerved data of OORT and 
RASMUSSEN (1971) up to 20 km. Values of V based on observations are almost 
nonexistent above about 20 km. The observed v values were thus extrapolated 
to 60 km by the use of simple analytic functions that matched the 20-km values 
and could be adjusted until the mean circulation resembled that of other 
Investigators, in particular that from the three-dimensional model of 
CUNNOLD, et al. (1975). 

The eddy fluxes are parameterized by the eddy diffusion coefficients 

K , K , and K . (The local "y-axis" is tangent to the meridian, i.e., 
vy yz zz 

dy ■ Rde). The K values were obtained by trial and error until reasonable 
agreement was obtained between calculated and observed distributions of 
Carbon 14 (JOHNSTON e£ al. , 1976) and the meridional distributions of the 
ozone column collated from reports of radiosonde observations by WILCOX et al. 
(1977). Except as otherwise noted, the eddy coefficients are those described 
in WHITTEN et al. (1977). 

5. Comparisons of the prediction and observations of the ozone distributions 

Figures 3 and 4 show the superposition of the BUV data (shown with ilOZ 
error bars) and the model predictions for the four seasons. The dramatic 
changes in the latitude dependence of the data with pressure-altitude during 
a single season are evident. At pressure-altitudes of 1 and 2 mb, there is 
also a dramatic change of the latitude dependence with season. However, at 
5 and 10 mb, the latitude dependence of the data is nearly independent of 
season. The model predictions of the changes in the latitude dependencies 
of the ozone abundances with season and with altitude agree well with the 
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observations. Although the agreement is good (note that linear rather than 
log scalea are used in Figs, 3 and 4), there are discrepancies thet are so 
much greater than the variability of the data that they require explanation. 
For example, the high-latitude predictions, particularly near the winter pole, 
are usually above the measurements. Although there is the possibility that 
the BUV experiment was somewhat less accurate at high latitudes because of 
the small amount of solar ultraviolet radiation to be bsckscattered and 
because of the obliqueness of the Incoming radiation, it seems more likely 
that the model parameters that are used in the predictions are less certain 
at these latitudes. 

To determine how the comparisons would be affected by other choices of 
atmospheric parameters that were assumed in the model or that were not well 
known, model runs were made for two temperatures and two choices of the 
mixing ratios of water vapor, odd chlorine and odd nitrogen. Two runs were 
also made to study the effects of changing the assumed transport parameters. 
The results of the parametric study are presented in the following sections. 

Sensitivity cf the comparisons to the assumed temperature fields 

To determine the sensitivity of the comparison to the temperature 
uncertainties, the model temperatures were uniformly Increased by 10*C at 
altitudes above 20 km and a new set of ozone predictions were generated. 

These predictions are plotted as dashed curves in Fig. 5, From this figure 
it can be seen that the predictions are lowered from 5X to 15X, depending on 
pressure height and latitude. In general, the Increased temperatures cause 
the predictions to fall below the observations and lessen the agreement 
between the observations and predictions. 
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The temperature fields that are used In the aodel for altitudes above 
20 km are derived from the rocket-data collation of NASTROM and BELMONT (1975); 
they represent long-term averages for each aaaaon. It Is unlikely that 
seasonal-average temperatures that occurred during the 1969-1971 period 
differed by as much as 10*C from thoae used in our model. In fact, the 
temperature measurements made by the Selective Chopper Radiometer, which were 
reported by BARNETT (1974) for the period from November 1970 to November 1971 
for pressure heights between 2 and 5 mb, are within ±3*C of those we used. 

If we assume that the temperatures in the stratosphere did not vary from those 
used in the model by more than ±3°C, then the effect of such a temperature 
uncertanlty is less than ±5Z at all latitudes and altitudes. 

Sensitivity of the comparisons to the assumed water vapor abundance 

The predicted ozone concentrations are sensitive to the amount of water 
vapor present because the water vapor reacts with 0(*D) atoms to produce odd 
hydrogen species which catalytically destroy ozone at high altitudes. Unfor- 
tunately, tonally averaged water-vapor measurements are not available as a 
function of latitude and season for the 1969-1971 period of the BUV observa- 
tions. Therefore, two different model assumptions about the amount of water 
vapor present were made to determine the sensitivity of the comparison to 
the range of water vapor likely to have been present during the observation 
period. 

A collation of stratospheric measurements of HjO at midlatitudes by 
HARRIES (1976) indicated that the average H 2 O mixing ratio varied from 3 ppmv 
at 15 km to 5 ppmv at 45 km with an rms deviation of 1.5 ppvm. The H 2 0 mea- 
surements of HILSENRATH £t al. (1977) show that the stratospheric water abundance 
is quite lnhomogenous in that water-vapor mixing ratios as low as 1.6 ppmv and 
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as high as 13 ppmv were observed In their aircraft surveys between 70* N and 
50* S latitude. 

When the assumption of the veter-vapor abundance In the model la changed 
from 4 ppmv to 8 ppmv, the model predicts a latitude-dependent osone decrease 
that is most pronounced at the upper altitudes (See Fig, 6 ). The effect 
is greatest at low latitudes and varies from 3% at 5 mb to 15% at 1 mb pressure 
altitude. 

Sensitivity of the comparisons to the assumed odd ehlonne abundance 

The effect on the predicted ozone profiles of the amount of odd chlorine 
(odd chlorine = Cl + HC1 + CIO + CINO 3 ) present is shown in Fig. 7, The 
solid curve shows the effects of assuming an asymptotic mixing ratio of 
1.6 ppbv, where the dashed curve shows the effects of assuming 3.2 ppbv of 
odd chlorine. When the lower value of the odd chlorine mixing ratio is 
used, the model predictions for HC1 are in agreement with the 15-to-20-km- 
altitude measurements of FARMER e_t al . (1976), ACKERMAN et «1. (1976), 

WILLIAMS et al . (1976), and EYRE and ROSCOE (1977). However, the predicted 
HC1 curves are a factor of two to three below the measurements of HCl at 
25 km and above. Doubling the assumed odd-chlorine mixing ratio to 3.2 ppbv 
in the model brings the model predictions of HCl into better agreement with 
the HCl measurements above 25 km. (See Fig. 8 .) However, an assumption of 
an asymptotic mixing ratio of at least 4 ppbv odd chlorine is needed to 
match the Cl and CIO measurements of ANDERSON £t al. (1977) and MENZIES 
(1979). Although a model assumption of 3.2 ppbv of odd chlorine does not 
significantly change the overall agreement between the model predictions 
and the measurements of ozone, it does result in better agreement with the 
measurements of HCl above 20 km and with Cl and CIO. 
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Sensitivity of the comparisons to the assumed odd nitrogen abundance 

To determine the sensitivity of the ozone comperlaon to the amount of 
odd nitrogen aaaumed to be preaent In the stratosphere, the model vaa run for 
two different aaaumptlona about the background level of odd nitrogen. The 
model calculates 13 ppbv of odd nitrogen at 40-km altitude when simulating 
the ambient atmosphere. This result Is nearly Independent of the latitude 
but It does depend on the altitude. For the aenaltlvity teat considered here, 
the odd-nltrogen mixing ratios at all latitudes and altitudes were multiplied 
by 1.5 and then the model was run for one simulated year. The 13 ppbv and 
the 18 ppbv of odd nitrogen used in this test are above the 11 ppbv measured 
by EVANS et al . (1976) and the 13 ppbv deduced by ACKERMAN et al . (1975) , 
but they are below the odd-nltrogen abundance implied by the NO 2 measurement 
of HARRIES jit al. (1976). The predicted effects of the Increased odd nitrogen 
are shown In Fig. 9. An examination of this figure shows that the predicted 
effects are negligible at 1 mb and small at 2 mb, but substantial at 5 and 
10 mb. That the predicted effects are more pronounced at lower altitudes than 
at higher altitudes shows that the destruction of ozone Is dominated by odd- 
nltrogen reactions at the lower altitudes and by other (l.e., odd hydrogen) 
reactions at the upper altitudes. It Is clear from the figure that changes 
In the model assisaptlons over this range of values can cause a 10Z change In 
the predicted ozone abundance at 5 mb. 

Sensitivity of the comparisons to the transport parameterization 

Because the comparisons are at altitudes where ozone is close to photo- 
chemical equilibrium, the co^iarlsona are not expected to be strongly Influ- 
enced by the exact choice of the transport parameters. The dashed curves In 
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Fig. 10 About comparisons for ths csss vhcrs both the winds and sddy transport 
wars reduced by a factor of two. This reduces the osona fluxes by a factor 
of two becauae neither the osona abundances nor thalr gradients are appreciably 
changed. From thla figure It can be seen that only alnlaial changes of the 
predictions at the 1- and 2-mb pressure heights occur except near the winter 
pole. However, significant Improvements in the comparisons at the 5- and KHsb 
heights are evident at both mid and high latitudes. In the polar regions, the 
agreement between the predictions and the observations Is significantly 
enhanced by the reduction of the transport parameters. The substantial changes 
that occur at all altitudes near the winter pole show the lmportancs of the 
transport parameterization when the ozone that is present is primarily that 
which has been transported from lower latitudes. In general, however, we can 
conclude that the comparisons at low and midlatitudes are, as expected, not 
highly sensitive to simple scaling changes in the transport parameters. 

6. Sunnary 

Model predictions of the variation of the ozone distribution with altitude, 
latitude, and season represent the state of our understanding of the mechanisms 
that control the production, loss, and movement of ozone. However, there are 
many uncertainties concerning the exact values of boundary conditions, input 
parameters, and the quantitative representation of transport, especially 
above balloon attitudes. Hence, the validation of atmospheric models Is 
critical to developing the capability of predicting the impact of anthro- 
pomorphic activities. The present comparison between model predictions and 
BUV ozone observations has shown that the predicted variation of the latitude 
dependence with altitude and with season agrees well with the observed 
variations. 
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The sensitivity of the model predictions to uncertslnites In the strato- 
spheric temperatures, the abundance of water vapor, the abundance of odd 
chlorine, the abundance of odd nitrogen, and the transport parameter laatlon 
was examined. The results of the sensitivity test Indicate that: 

1. A ±3°C temperature uncertainty contributes a i5Z change in the 
predicted ozone abundances. 

2. When the assumed stratospheric vater-vapor mixing ratio is changed 
from 4 ppmv to 8 ppmv, the model predictions change from 32 to 152, depending 
on both latitude and altitude. 

3. A model assumption of an asymptotic mixing ratio of 3.2 ppbv for 
odd chlorine lowers the predicted ozone abundance by 152 but does not 
significantly chauge the overall agreement with the observed ozone abundances. 

4. A 502 increase in the assumed odd-nitrogen level causes the pre- 
dictions of the ozone abundance to decrease by 1C2 at 5 mb. 

5. Although the model predictions of the ozone abundances at these 
altitudes are not highly sensitive to the exact values of the transport para- 
meters, the agreement between the predictions and the observations is 
significantly improved at mid and high latitudes when the transport parameters 
are reduced by a factor of two. 
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APPENDIX 


Photochemical Reaction Rataa 



Chemical react Iona and corresponding rate 

coefficients 

No 

• Reaction 

Rate coe. ficient* 

Reference 

1 

0 ♦ 0 3 -*• 0 2 ♦ o 2 

1.5*10" r exp(-2218/T) 

b 

2 

o + o 2 + m-*o 3 + m 

1.1*10“ j4 exp(510/T) 

c 

3 

no 2 + 0 NO -f 0 2 

9.3*1C"12 

b 

4 

NO + 0 3 -*• N0 2 ♦ 0 2 

2. 3* 10“ 12 exp(-1450/T) 

b 

5 

N0 3 + NO -*> N0 2 + N0 2 

8.7<10“ 12 

c 

6 

0 3 + N0 2 N0 3 + 0 2 

1. 2x10“ 13 exp(-2450/T) 

d 

7 

0( ! D) +M+0+M 

1. 659*10“* 1 exp(107/T) 
+ 0.58 10“* 1 exp(67/T) 

d 

8 

0 + OH -► H + 0 2 

4.0*10 -11 

b 

0 

0 3 + H -*■ OH + 0 2 

1.4*10“ 10 exp(-470/T) 

b 

10 

NO + 0 + M -*• N0 2 + M 

1.55*10“ 32 exp(584/T) 

d 

11 

N0 2 + N0 3 + M ■+ N 2 0 r + M 

e 

TURC0 and WHITTEN (1975) 

12 

0 ( *D) + H 2 0 20H 

2. 3* 10“ 10 

d 

13 

o 3 + ho 2 -► oh + :j 2 

1.1*10“ 14 exp(-580/T) 

ZAHNISER and HOWARD (1978) 

14 

0 + H0 2 •+ OH + ' 2 

3. 5*10“ 11 

d 

15 

N0 2 ♦ OH + M -*• HN0 3 + M 

f 

d 

16 

OH + HN0 3 H 2 0 + N0 3 

8. 5* 10- 14 

b 

17 

0 + 0 + M * J 2 + M 

3*10’ 33 (300/T) 3 

c 

18 

OH + H0 2 - H 2 0 + 0 2 

4x10“ 11 

b 

19 

OH + CO - H ♦ C0 2 

1. 35*10“ 1 3 (1+P ) 

atm 

b 

?0 

H + 0 2 - M -* H0 2 + M 

2. 1*10“ 32 exp(290/T) 

d 

See 

footnc s at end of table, 

D. 
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APPENDIX - Continued 


No. 

Reaction 

Rate coefficient* 


Reference 


21 

0 + HN0 3 ■* OH + N0 3 

1*10“ 1,1 

c 




22 

OH + 0 3 -*• H0 2 + 0 2 

1.6*10" 12 e*pl(-940/T) 

b 




23 

NO + H0 2 ■* N0 2 + OH 

3.4*10" 12 expl(250/T) 

b 




24 

N0 2 + 0 ♦ M N0 3 + M 

g 

TURCO 

and 

WHITTEN 

(1975) 

25 

N 2 0 5 + M -♦ N0 2 + N0 3 + M h 

TURCO 

and 

WHITTEN 

(1975) 

26 

H + OH + M -*> H 2 0 + M 

6.1*lfT 26 T -2 

c 




27 

N + 0 3 -*■ NO + 0 2 

5*10~ 12 exp(-650/T) 

d 




28 

N + NO -*■ N 2 + 0 

3. 4* 10" 1 1 

b 




29 

N + 0 2 -*■ NO + 0 

5. 5* 10" 12 exp(-3220/T) 

d 




30 

N + OH -*• NO + H 

5. 3*10" 11 

c 




31 

NO + OH + M •* HN0 2 + M 

2 . 2* 10" 3 2 «xp(1100/T) 

c 




32 

OH + OH H 2 0 + 0 

1x10“ 11 exp(-500/T) 

b 




33 

OH + HN0 2 -♦ H 2 0 + N0 2 

3*10" 14 

c 




34 

0( ] D) + N 2 + M -*• N 2 0 ♦ 

M 3.5*10" 37 

d 




35 

0( J D) + N 2 0 N 2 + 0 2 

5. 1*10" 11 

b 




36 

0{'D) + N 2 0 -*• NO + NO 

5.9*10" u 

b 




37 

0 HN0 2 -*■ OH + N0 2 

1.5*10" 14 

c 




38 

OH + OH + M H 2 0 2 + M 

1.25*10" 32 exp(900/T) 

d 




39 

OH + H 2 0 2 -*• H 2 0 + H0 2 

1*10" 11 exp(-750/T) 

d 




40 

ho 2 + H0 2 H 2 0 2 + 0 2 

2.5*10" 12 

d 




41 

0 + H 2 0 2 -*■ OH 4 H0 2 

2.8*10" 12 «xp(-2125/T) 

b 





See footnotes at end of table, p. 24 . 
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APPENDIX - Continued 


No. 

Reaction 

Rate coefficient* 

Ref trance 

42 

HC1 ♦ OH H 2 0 ♦ Cl 

2. 8x10“ 12 exp(-425/T) 

ZAHNISER at al. (1974) 




SMITH and ZELLER (1974) 




RAV1SHANKARA at al. (1977 

43 

HC1 + 0 OH + Cl 

1.14x10” 11 exp(-3370/T) 

b 

44 

Cl + CHi* ♦ HC1 ♦ CH 3 

9. 9x10“ 12 exp(-1359/T) 

b 

43 

Cl + H 2 •* HC1 + H 

3.5xl0" u txp(-2290/T) 

d 

46 

Cl + 0 3 CIO + 0 2 

2. 8x10” 11 exp(-257/T) 

b 

47 

CIO ♦ 0 Cl + 0 2 

7. 7x10’ 11 exp(-130/T) 

d 

48 

CIO + NO -*• Cl + N0 2 

7.8xl0* 12 exp(250/T) 

b 

49 

CIO + CO Cl + CO 

1.7xlO"15 

c 

50 

CIO + N0 2 +M-* C10N0 2 + M 

1 

d 

51 

C£0N0 2 + 0 ♦ CIO ♦ N0 3 

3x10’ 12 exp(-808/T) 

b 

52 

Cl + H0 2 ■* HCl ♦ 0 2 

4.5xl0“ n 

b 

53 

0( J D) + H 2 ♦ OH + H 

9.9*10-11 

d 

54 

O^D) + CH 4 -*• OH ♦ CH 3 

1.3x10”! 0 

d 

55 

OH + CH 4 -♦ H 2 0 ♦ CH 3 

2. 4x10' l 2 exp(-1710/T) 

b 

56 

CH 3 + 0 2 +M-»CH 3 0 2 + M 

J 

TURCO and WHITTEN (1975) 

57 

CH 3 + 0 - CH 2 0 + H 

1. 2“10 _1 0 

k 

58 

2CH 3 0 2 -*■ 2CH 3 0 + 0 2 

2. 0x10’ 1 5 

k 

59 

CH 3 0 2 + NO CH 3 0 ♦ N0 2 

8x10“ 12 

b 

60 

ch 3 o 2 ♦ ho 2 *ch 4 o 2 ♦ o 2 

lxlO” 12 

b 

61 

ch 4 o 2 ♦ oh-*ch 3 o 2 ♦ h 2 o 

6. 2xl0-l 2 exp(-750/T) 

b 

See 

/ootnotee at end of table 

, P- 24 - 
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No. 

Reaction 

Bate coefficient* 


Reference 

62 

ch 3 o * o 2 ■* ch 2 o ♦ ho 2 

5*10" 13 exp(-2000/T) 

b 


63 

CH 2 0 ♦ 0 -► CHO ♦ OH 

2.8*10" n «xp(-1540/T) 

b 


64 

CH 2 0 + OH -► CHO ♦ H 2 0 

1.7*10" n exp(-100/T) 

b 


65 

CHO + 0 2 ■» CO + H0 2 

5*10" 12 

b 


66 

CHO + 0 CO + OH 

1.0«10" 10 

k 


67 

CHO + 0 *- C0 2 + H 

7.0*10 -11 

k 


68 

CH3 CL+ oh-h 2 o + ch 2 o 

2. MO" 12 exp (-1142 /T) 

b 



+ CLO 




69 

CH 2 0 + CL HCL + CHO 

7.5*10“ n 

b 



Photodlssoclatlon Rates 


No. 


Reaction 


Dissociation 

rate* 


Reference* 


Jl 0 2 ♦ hv - 20( 3 P) 9.7M10" 6 


J2 

0 3 ♦ hv 0 2 4 0( 3 P) 

3.47*10"** 

J 3 

0 3 ♦ hv -* 0 ? *6 0( 1 D) 

5.06*10" 3 

J*. 

NO + hv -* N ♦ 0 

2.78*10"* 

Js 

N0 2 + hv •* NO 4 0 

5.88*10" 3 


See footnote* at end of table, p. 24. 


HUDSON and HAHLE (1972); OGAVA 
(1968); BLAKE at al. (1966); 
VATANABE «t al. (1953) 

INN and TANAKA (1953); CRIGGS 
(1968); JONES and WAYNE (1969) 
INN and TANAKA (1953); CR1GCS 
(1968); JONES and WAYNE (1969) 
CIESLIK and NICOLET (1973) 

HALL and BLACET (1952); NAKAYAHA 
et al. (1959); PITTS et al. 
(1964) 
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APPENDIX - Continued 


No. 

Reaction 

Dissociation 

rata 1 

Raferancas 

H 

H 2 0 ♦ hv •» OH + H 

1 . 6 * 10 “® 

WATANABE and ZELIKOFP (1953) 

J 7 

HN0 3 + hv ♦ N0 2 ♦ OH 

7. 17*10“ 5 

JOHNSTON and CRAHAM (1974) 

Ja 

HN0 2 ♦ hv NO ♦ OH 

4.29*10“** 

JOHNSTON and CRAHAM (1972) 

J 9 

N 2 0 ♦ hv -► N 2 4- 0(*D) 

4.54*10“ 7 

JOHNSTON and SELVYN (1975) 

J 10 

N0 3 + hv -* N0 2 ♦ 0-*N0 ♦ 

0 2 1 . 66 * 10” 2 

JOHNSTON and GRAHAM (1974) 



2 . 66 * 10“ 2 


J 11 

N 2 0 5 ♦ hv ■* 2N0 2 + 0 

3.84*10“** 

GRAHAM (1975) 

J 12 

H0 2 + hv -*• OH + 0 

5.43*10-** 

PADKERT and JOHNSTON (1973) 

J 1 3 

H 2 0 2 ♦ hv 20H 

7.64*10“ 5 

MOLINA at al. (197?) 

Jl4 

HC£ ♦ hv -*■ H + Cl 

5.65*10“ 6 

MYER and SAMSON (1970) 

J 15 

CIO ♦ hv Cl + 0 

3.36*10“ 3 

JOHNSTON at al. (1969); 

» 




LANGHOFF at al. (1977) 

J lt 

C!0N0 2 ♦ hv -► CIO ♦ N0 2 

3.52*10“** 

b 

J 17 

CF 2 C1 2 + hv -► 2C1 

9.92*10“ 7 

■ 

J 18 

CFCtj ♦ hv -*• 2. SCI 

6.27*10“* 

■ 

J 19 

CH 2 0 ♦ hv -* CHO ♦ H 

5. 0*10“ 5 

TURCO (1975) 

J 20 

CH 2 0 + hv CO + H 2 

7.0*10“ 3 

TURCO (1975) 

J 21 

CH,,0 : ♦ hv -► CH 3 0 ♦ OH 




CH 3 CI ♦ hv -*• CH 3 ♦ Cl 

1.05*10“ 7 

ROBBINS (1976) 


*In units of ca 3 see -1 for blaaoleculsr reset Ions and ca* sac -1 for 


trlaoleculsr reactions. 

b NASA Panel for Data Evaluation, 1979. 
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APPENDIX - Continued 

"» - ■ j imu-j-ag— — — a— .11, mjl ■■;-TTr 

c Rate coefficients are those recoaBcnded In HAMPSOM and CAXVIN (1975), 

^Rate coefficients are the valuea recoaBcnded In HUDSON (1777;. See 
Table 1 in HUDSON (lfp7) and corresponding notes therein for a discussion of 
the pertinent laboratory aeasuraaents upon which these values are based. 

*The three-body reaction-rate coefficient t 
ki i • (1.9*10" 30 )/(1 4 5V ) • (3 4 5V )/(3 4 V,.) ca* sec" 1 , where V - 5»10* l 9 [M]. 

f The equivalent three-body reaction-rate coefficient for the formation 
of HNO 3 is ti'’en by: kjj ■ (0.94/M) (280/T) 1 ^ 2 10 C a 6 sec" 1 , where: 

A • 31.62273 - 0.258304 • 2 - 0.0889287 • Z 2 4 2.520173*10“ 3 • Z 3 
B • -327.372 4 44.5586 • Z - 1.38092 • Z 2 
Z • logjolO.78 • M] 

T ■ teaperature 

*The reaction-rate coefficient k 2 t, is given by 
k 2 4 - 2.8*10" 31 /(1 4 2.55-10 " 20 M) cm 6 sec" 1 . 

^The reaction-rate coefficient for the colllslonal destruction of N 2 0 ^ is 
given by k 25 - 2.2*10 " 5 V q (4 »• Vj)/I(1 4 Vj)(4 4 Vj/ 6 )), where V fl - exp(-9700/T) 

and Vj - 2.3*10 " 19 exp(900/T) M. 

i © 

The foraatlon rate of chlorine nitrate is givtn by: 

k^o “ 3. 3*10 ” 2 3 t" 3 • 3 *• / < 1 4 8.7-10 " 9 T " 0 * 6 M 0 * 5 ) c*. c sec" 1 . 

^The pressure-dependent three-body reaction-rate coefficient is 

k Sfe - 2.6*10‘ 31 /(1 4 6.047-10" 19 • M) . 

k See POPPOFF tt si. (1978) for reference list and detailed dlscuasion of 

the aethane oxidation cycle. 

^Photodlseoclation rates represent dlurnelly sveraged vslucs (altitude 60 ka, 
latitude 45* N, Autumn). The units are sec* 1 . 
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APPENDIX - Concluded 


Reaction rate* are calculated front the cross sections tabulated in 
HUDSON (1977). 

n Slnce no photodlssoclatlon cross sections are available at the time of 
writing, we have assumed (after TURCO and WHITTEN, 1975) that the dissociation 
rate of CH4O2 Is comparable to that of H2O2. 
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Figure Captions 


1. Schematic diagram of the BUV experiment. 

2. Effective scattering levels for solar radiation in the nadir direction 
of the satellite for all orders of scattering; solar zenith angle ■ 60* 
and total ozone ■ 336 m atm- era. (From HEATH et^ jal. , 1973.) 

3. Comparison of the Goddard BUV observations with the Ames two-dimensional- 
model predictions for fall and winter. The symbols represent the BUV 
data for 1970 and 1971 while the solid curve represents the model 
predictions. 

4. Comparisons of the Goddard BUV observations with the Ames two-dimensional- 
model predictions for spring and summer. The symbols represent the BUV 
data for 1970 and 1971 while the solid curve represents the model 
predictions. 

5. Comparison of the model predictions and BUV data where the model 
temperatures are Increased by 10*C. 

6. Comparison of the model predictions and the BUV data when the water vapor 
mixing ratio in the model is changed from 4 ppmv to 8 ppmv. 

7. Comparison of the model predictions and the BUV data when the odd-chlorine 
mixing ratio in the model is changed from 1.6 ppbv to 3.2 ppbv. 

8. Comparison of the predicted and measured HC1 abundance in the fall when 
an odd-chlorine asymptotic mixing ratio of 3.2 ppbv is assumed. 

9. Comparison of the model predictions and the BUV data when the odd- 
nilrogen mixing ratio in the model is multiplied by 1.5. 

10. Comparison of the model predictions and the BUV data when the model 
transport parameters are divided by two. 
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